Auxiliary Variable-based Balancing (AVB) for source term treatment in open channel simulations

6 Practical engineering applications of open channel flow modelling involve geometric terms 7 arising from variations in channel shape, bottom slope and friction. This paper presents the 8 family of schemes that satisfy the generalised C-property for which static equilibrium is a 9 particular case, in the framework of one-dimensional open channel flows. This approach, 10 named Auxiliary Variable-based Balancing, consists of using an auxiliary variable in place of 11 the flow variables in the diffusive part of the flux estimate. The auxiliary variable is defined 12 so as to achieve a zero gradient under steady-state conditions, whatever the geometry. Many 13 approaches presented in the litterature can be viewed as a particular AVB case. Three auxil14 iary variables are presented in this paper: water elevation, specific force and hydraulic head. 15 The methodology is applied to three classical Riemann solvers: HLL, Roe and the Q-scheme. 16 The results are compared on five test-cases: three steady-state configurations including fric17 tion, singular head losses and variations in bottom elevation, channel width and banks slope 18 and two transient test-case (dam-break problems on rectangular and triangular channel). 19 In each case, the auxiliary variable that best preserves the steady-state configuration is the 20 hydraulic head. Besides, using the head as auxiliary variable allows head loss functions due 21 to singularities to be incorporated directly in the governing equations, without the need 22 for internal boundaries. However, it is generally less accurate when sharp transients are 23 involved. 24